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ABSTRACT 

Context. Precise determinations of the chemical composition in early B-type stars constitute fundamental observational constraints on 
stellar and galactochemical evolution. Carbon, in particular, is one of the most abundant metals in the Universe but analyses in early- 
type stars are known to show inconclusive results. Large discrepancies between analyses of different lines in C ii, a failure to establish 
the C ii/iii ionization balance and the derivation of systematically lower abundances than from other indicators like H ii regions and 
young FG-type stars pose long-standing problems. 

Aims. We discuss improvements to the non-LTE modelling of the visual line spectrum and to the spectral analysis of early B-type stars, 
as well as their consequences for stellar parameter and abundance derivations. The most relevant sources of systematic uncertainies and 
their effects on the analysis are investigated. Consequences for the present-day carbon abundance in the solar vicinity are discussed. 
Methods. We present a comprehensive and robust C ii/iii/iv model for non-LTE line-formation calculations based on carefully selected 
atomic data. The model is calibrated with high-S/N spectra of six apparently slow-rotating early B-type dwarfs and giants, which cover 
a wide parameter range and are randomly distributed in the solar neighbourhood. A self-consistent quantitative spectrum analysis is 
performed using an extensive iteration scheme to determine stellar atmospheric parameters and to select the appropriate atomic data 
used for the derivation of chemical abundances. 

Results. We establish the carbon ionization balance for all sample stars based on a unique set of input atomic data. Consistency 
is achieved for all modelled carbon lines of the sample stars. Highly accurate atmospheric parameters and a homogeneous carbon 
abundance of log(C/H)-l- 12 = 8.32 ±0.04 with reduced systematic errors are derived. Present evolution models for massive stars 
indicate that this value may require only a small adjustment because of the effects of rotational mixing, by < -1-0.05 dex per sample 
star. This results in a present-day stellar carbon abundance in the solar neighbourhood, which is in good agreement with recent 
determinations of the solar value and with the gas-phase abundance of the Orion H ii region. Our finding of a homogeneous present- 
day carbon abundance also conforms with predictions of chemical-evolution models for the Galaxy. Moreover, the present approach 
allows us to constrain the effects of systematic errors on fundamental parameters and abundances. This suggests that most of the 
difficulties found in previous work may be related to large systematic effects in the atmospheric parameter determination and/or 
inaccuracies in the atomic data. 

Key words. Line: formation - Stars: early type, fundamental parameters, abundances - Galaxy: abundances, solar neighbourhood 



1. Introduction 

Studies of the chemical composition of early-type stars provide 
valuable observational constraints to our understanding of broad 
fields like stellar evolution, nucleosynthesis and galactochemical 
evolution. The accuracy of elemental abundance determinations 
defines the degree to which theory can be tested. However, quan- 
titative spectroscopy relies on many model assumptions itself. A 
detailed understanding of the interaction between the radiation 
field and the plasma in the stellar atmosphere is thus required. 
Any weakness in the physical model limits the analysis even 
when high quality in observation and data reduction is achieved. 

Concerning Galactochemical evolution (e.g. models of Hou 
et al. 120001 Chiappini et al. 120011 [2003 ), early-type stars are 
important probes for the current chemical composition of the 
Galaxy. They provide spatial information of elemental abun- 
dances also beyond the solar vicinity, allowing Galactic metal- 
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licity gradients to be derived (e.g. Gummersbach et al. '19981 
Rolleston et al. l2000l Daflon & Cunha l2004b . Abundances from 
young stars complement data from H n regions, allowing the cur- 
rent status of the nucleosynthesis of metals in the course of the 
cosmic cycle of matter to be constrained. Pristine abundances 
can be derived not only for Galactic early-type stars, but also 
for objects in the metal-poor environments of the Magellanic 
Clouds, using the present generation of large telescopes and 
high-resolution spectrographs (e.g. Korn et al. 2002, 2005) 
Rolleston et al. 120031 Hunter et al. I2007l i. Such investigations 
can be extended to more distant galaxies of the Local Group, 
and even beyond, by analysing intermediate-resolution spectra 
of early-type supergiants (e.g. Trundle et al. 120021 Urbaneja et 
al. 2005ab). Moreover, CNO abundances derived from quanti- 
tative spectral analyses of massive early-type stars provide fun- 
damental observational constraints for stellar evolution models. 
Basic stellar parameters and abundance patterns of the light el- 
ements resulting from mixing with nuclear-processed matter fa- 
cilitates an empirical evaluation of different evolution models 
(e.g. Heger & Langer |2000l Maeder & Mevnet |2000] l. 
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Carbon plays a special role among the light elements. It is 
a primary element created in the fundamental Sa-process and as 
such it provides the seed for the subsequent synthesis of all heav- 
ier elements (Burbidge et al. |1957t Cameron |1957l l. Carbon is an 
essential catalyst for the nucleosynthesis of H into He through 
the CNO cycle in massive and intermediate-mass stars. The ele- 
ment also constitutes the basis of all organic chemistry. 

Carbon abundances derived from early-type stars were the 
subject of numerous spectroscopic studies in the past decades, 
from the pioneering work of Unsold (1942) to numerous others 
since then. The quality of observed data and the complexity of 
the model calculations and spectral analysis have been improved 
with great effort over these years, however some inconsistencies 
remained until the present. The 'classical' carbon problem in 
early-type stars addresses the extremely low abundances derived 
from the strong C ii /14267 A doublet when compared with those 
derived from weaker lines, as pointed out by Kane et al. ( 1980), 
referring to work by Hardorp & Scholz ( 1970) and Kodaira & 
Scholz d 19701 1. Later on, more complications were found, such 
as discrepancies from both strong C ii /1/16578/82 A and 4267 A 
multiplets, /1/16578/82 A and the weak lines, among the C ii weak 
lines themselves and in abundances derived from C ii and C iii. 

In nearby, bright and apparently slow-rotating stars it is 
possible to derive carbon abundances from weak lines which 
are less sensitive to the details of the analysis. However, stud- 
ies of extragalactic or fast-rotating stars are restricted to the 
stronger lines because of S/N-constraints or rotational smearing. 
In those cases, one might be interested to analyse not only the 
C n /1/16578/82 A and the 4267 A multiplets, but also other strong 
C m or C IV Unes, in order to assure ionization equilibrium. 

A crucial step in the spectral modelling was to abandon 
the approximation of local thermodynamic equilibrium (LTE) 
in line-formation calculations and allow for deviations (non- 
LTE). Several model atoms have been discussed in literature 
(Lennon 1983; SakhibulHn 1987; Eber & Butler [19881 Grigsby 
et al. [1992; Sigut 1996; Przybilla et al. IIOOTI Lanz & Hu- 
benv l2003. .2007 ). The model atom of Eber & Butler found wide 
appUcation for abundance analyses of mostly unevolved early- 
type stars in the solar neighbourhood (e.g. Gies & Lambert[1992j 
KiUan [T992l Cun ha & L ambert 199 4; Gum mersbachet al. [T998l 
Andrievsky et al. [19991 Daflon et al. [19991 2001ab). 

The C 11/14267 A multiplet poses a challenge to non-LTE 
line-formation calculations (e.g. Lambert 119931 Sigut 1996). 
Even in non-LTE studies this and the /1/16578/82A multiplet 
failed to reproduce observation consistently (Grigsby et al. II 9921 
Hunter et al. |2007l l. Systematically lower abundances are derived 
than from other weaker C ii lines (Gies & Lambert 1992), as in 
the LTE case. Observed trends for C iii lines may also be poorly 
matched by non-LTE model calculations (Grigsby et al. 119921 ). A 
failure to establish the C ii/iii ionization equilibrium in non-LTE 
was also reported. Differences in abundance from the two ions 
can amount up to a factor ~5-10 (Daflon et al. 2001b; Hunter 
et al. 12007 ). There are even more inconsistencies with published 
carbon abundances from early-type stars in a broader context 
that require other explanation, as most of the published studies 
avoid the lines sensitive to non-LTE effects. 

A comparison of the available studies of early B-type stars 
indicates that present-day carbon abundances in the solar vicin- 
ity are highly inhomogeneous (even for stars within a single 
cluster) and largely sub-solar This is in contrast to the find- 
ings of a uniform abundance in the gas-phase of the interstellar 
medium (ISM) within 1.5 kpc of the Sun (Sofia & Meyer [SOOTl 
and references therein). It cannot be understood from current 



stellar and Galactochemical evolution models either Young 
massive stars that form out of a molecular cloud within a short 
timescale can be expected to show a homogeneous carbon abun- 
dancfl which also coincides with that of their surrounding H ii 
nebula. Moreover, this abundance is expected to be higher than 
that of objects from previous generations of star formation in 
their neighbourhood, like young (<2 Gyr) F and G stars, and that 
of the Sun. In reality, significant systematic differences exist (see 
e.g. Sofia & Meyer jlOUTl Herrero [2003l l. shedding doubt on the 
reliability of carbon abundances derived from early B-type stars. 

Quantitative spectroscopic studies depend on the quality of 
the observational data, the physical approximations made in the 
modelling approach (stellar atmospheres, model atoms) and the 
analysis technique. The analysis procedure also involves several 
steps (data reduction, stellar parameter determination and abun- 
dance derivation), which are all critical for the quality of the final 
results. It is difficult to control systematic effects in the whole 
process. Therefore, the sources of discrepant results from the lit- 
erature are hard to be identified. 

One underestimated problem in the derivation of the chemi- 
cal composition is certainly the selection of the input atomic data 
for non-LTE line-formation calculations (see e.g. Sigut[T996|for 
C ii). Even when the model atmosphere computations reproduce 
the plasma conditions in the stellar photosphere realistically, use 
of different atomic data for the line-formation calculations can 
result in large systematic errors in the final abundances. The 
structure of atoms and therefore the atomic line spectra are gov- 
erned by the laws of quantum mechanics, i.e. they are essentially 
independent of the environment the atoms are located in. This 
property can be used to calibrate model atoms for non-LTE cal- 
culations, by demanding that the model should reproduce obser- 
vation reliably over the entire parameter space of relevance. In 
practice, this means that model atoms need to be tested for a few 
cases of well-established plasma conditions. The atmospheres 
of early B-type stars far from the Eddington limit provide such 
testbeds, as they can be described well by rather simple physical 
models. A small sample of stars, spanning the relevant range in 
spectral type and luminosity, is sufficient for such tests. 

The present study aims to solve several inconsistencies con- 
cerning the derivation of carbon abundances in early B-type stars 
and to apply the new modelling and analysis technique in stars 
of the solar vicinity. It focuses on both i) the critical selection 
of appropriate atomic input data from ab-initio calculations and 
approximation formulae and ii) a self-consistent derivation of 
the atmospheric parameters of the sample stars from ionization 
equilibria. This detailed and reliable study allows us to derive 
abundances with unprecedented accuracy for a so far small but 
very well constrained sample of stars covering strategic points 
of the stellar parameter space and randomly distributed in OB 
associations and in the field of the solar neighbourhood. 

This study is part of our efforts to improve on the accu- 
racy of quantitative spectroscopy in early B-type stars (Nieva 
& Przybilla [20021 Paper 1; Nieva & Prz ybilla [20061 Paper 2) 
and related objects (Przybilla et al. '2 0051 2006a-c). The hybrid 
non-LTE approach (LTE model atmospheres -i- non-LTE line for- 
mation) was shown to be adequate for analyses of OB-type main 
sequence and giant stars in Paper 1, allowing highly consistent 
results to be derived. In particular, practically all measurable H 
and He lines from the B aimer jump to the near-IR (where avail- 
able) have been reproduced with a single model for each of the 
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present sample stars, meeting a crucial boundary condition for 
all further analyses. Paper 2 discusses a solution to the 'classi- 
cal' carbon problem, i.e. discrepant carbon abundances from C ii 
/14267 A and weaker lines, and also a solution to the notorious 
discrepancy of the prominent /i/l6578/82 A and the 4267 A mul- 
tiplets. There, consistency was achieved for up to 21 C ii features 
of our sample stars. Those results are part of the present work, 
which discusses the topic in a much broader context. Here, we 
present for the first time a quantitative agreement of abundances 
from C II and C iii lines for this kind of star We also study C iv 
lines and establish the C ii/iii/iv ionization equilibrium in a self- 
consistent way, never considered before in similar work. The 
linelist rises up to 40 features. In the present work we inves- 
tigate the systematic errors arising due to uncertainties in the 
atomic data and stellar parameters and quantify the difference 
between LTE and non-LTE analyses, for the first time in this ex- 
tent. Implications on effective temperature scales for abundance 
determitations are discussed. This extensive study provides us 
with reliable carbon abundances for a discussion in the context 
of galactochemical evolution. 

The paper is organised as follows: Section 2 describes the 
model calculations including a description of the carbon model. 
Section 3 describes the empirical model calibration via an ex- 
tensive and self-consistent iteration and the sensitivity of carbon 
line-formation calculations to atomic data and atmospheric 
parameter variations. Section 4 summarises the final results for 
parameters and carbon abundances of individual lines for the 
sample stars. In Sect. 5 we provide a comparison of our results 
to previous studies and the conclusions on the present-day 
carbon abundance in the solar neighbourhood are discussed in 
Sect. 6. A summary is given in Sect. 7. 

2. Model calculations 

The non-LTE line-formation computations for carbon follow the 
methodology discussed in our previous study of the hydrogen 
and helium spectra in early-type stars (Paper 1). A hybrid 
non-LTE approach is employed to solve the restricted non-LTE 
problem on the basis of prescribed LTE atmospheres. This tech- 
nique provides an efficient way to compute realistic synthetic 
spectra in all cases where the atmospheric structure is close to 
LTE. The computational efforts can thus be concentrated on 
robust non-LTE line-formation calculations. 

2.1. Model atmospheres and programs 

The model atmospheres are computed with the Atlas9 code 
(Kurucz 1993b ), which assumes plane -parallel geometry, chem- 
ical homogeneity and hydrostatic, radiative and local thermo- 
dynamic equilibrium (LTE). Line blanketing is realised via 
consideration of Opacity Distribution Functions (ODFs, Ku- 
rucz [19933- We adopt standard solar abundances (Grevesse & 
Sauval [r998l l in all computations, as in Paper 1. Helium abun- 
dance and microturbulence are adjusted according to the results 
from the quantitative analysis of the programme stars - requir- 
ing some iteration, as explained below. The model atmospheres 
are held fixed in the non-LTE calculations. 

Non-LTE level populations and model spectra are obtained 
with recent versions of Detail and Surface (Giddings 119811 
Butler & Giddings [19851 both updated by K. Butler). The 
coupled radiative transfer and statistical equilibrium equations 
are solved with Detail, employing an Accelerated Lambda 
Iteration scheme of Rybicki & Hummer ( I199II ). This allows even 
complex ions to be treated in a realistic way. Continuous opaci- 



ties due to hydrogen and helium (for actual abundances) are con- 
sidered in non-LTE, line blocking is accounted for in LTE via 
Kurucz' ODFs. Synthetic spectra are calculated with Surface, 
using refined line-broadening theories. Microturbulence is con- 
sistently accounted for in both steps with Detail and Surface. 

Non-LTE level populations for hydrogen and He i/ii are com- 
puted using recent model atoms by Przybilla & Butler (I2004I I 
and Przybilla ( |2005 |), respectively. The synthetic H and Hei/n 
lines give excellent agreement with observations, in the optical 
as well as in the near-IR (Paper 1). This is a prerequisite for 
modelling metal lines which overlap with the (broad) hydrogen 
or helium features. In particular, the C ii /i/l6578/82 A doublet is 
affected in the present case. 

2.2. The C ii-iv model atom 

We give a short summary of the input atomic data for the con- 
struction of the C ii/iii/iv model atom in this section. A more de- 
tailed description of our motivation for choosing these atomic 
data will be given in Section 3.2, as this turned out to be critical 
for the realistic modelling of the observed spectra. 

Cii. This model ion considers LS' -coupled terms up to 
principal quantum number n = 10 and angular momentum ( = 9 
(66 levels) explicitly in the non-LTE calculations, with all 
fine-structure sub-levels combined into one. Additional lev- 
els up to n = 14 are computed in LTE relative to the ground 
state of Cm. Level energies are adopted from Moore ( I1993I I, 
Sigut d 1 9961 1 and Quinet (il998 ). The doublet and quartet 
spin systems are treated simultaneously. Oscillator strengths 
(^/-values) from three sources are considered: fine-structure 
data from ab-initio computations using the multiconfiguration 
Hartree-Fock method in the Breit-Pauli approximation of Froese 
Fischer & Tachiev (120041 FFT04), data from application of 
the Breit-Pauli /?-matrix method (Nahar 2002a N02a) and re- 
sults obtained in the Opacity Project (OP) from the /?-matrix 
method assuming LS-coupling (Yan, Taylor & Seaton 1987}. 
Our primary source of ^/-values is FFT04, followed by OP 
and N02a for the remaining transitions. Intercombination tran- 
sitions are neglected because of very small oscillator strengths. 
Photoionizations cross-sections are adopted from the OP (Yan 
& Seaton [T987l l for levels up to n = 9 and i-3 with a correction 
of the threshold frequencies to observed values and from Nahar 
& Pradhan (119971 NP97) for the remainder Effective collision 
strengths for electron impact excitation among the lowest 16 LS- 
states are adopted from /?-matrix computations of Wilson, Bell 
& Hudson ( 120051 |20071 i. An empirical increase by a factor two 
was applied to the 3s^S-3p^P° data (see Sect. 13.2b . Collisional 
excitation for transitions without detailed data are treated using 
the Van Regemorter ( 1962) approximation - in the optically al- 
lowed case - and via the semi-empirical Allen (I1973I I formula 
in the optically forbidden case. Collision strengths Q varying 
between 0.01 (An > 4) to 100 (An - 0) are employed, as 
suggested by evaluation of the detailed data from ab-initio com- 
putations of Wilson et al. (2005, 2007). Collisional ionization 
rates are evaluated according to the Seaton ( |I9621 l approxima- 
tion. Threshold photoionization cross-sections are adopted from 
OP and NP97, allowing for an empirical correction of one order 
of magnitude higher for the 6f ^F° and 6g ^G levels - correspond- 
ing to the upper levels of the C ii AA6 151 and 6462 A transitions, 
respectively. 

C m. This model accounts for LS -terms up to n = 7 and ( = 1 
(70 levels) explicitly in the statistical equilibrium calculations. 
In a similar way to Cn, levels up to n-l4 are computed in 
LTE relative to the ground state of the next ionization stage. The 
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two spin systems (singlet and triplet) are treated simultaneously. 
Level energies are taken from the NIST databas^ and energies 
for 9 of the highest levels are adopted from N02a. Two sources 
are considered for oscillator strengths: N02a and additional val- 
ues from Eber (,1987 ). Intercombination transitions are also 
implemented when their values are non-negligible (/>10 
Photoionization cross-sections are taken from S. Nahar's web- 
pag43 Maxwellian-averaged collision strengths for electron im- 
pact excitation among the lowest 24 terms are adopted from the 
/?-matrix computations of Mitnik et al. ( I2003I I. Collisional exci- 
tation for the remaining transitions and collisional ionization are 
treated in analogy to the Cii ion, using appropriate g/- values 
and threshold photoionization cross-sections. 

C IV. LS -terms up to n = 10 and £ = 9 (53 levels) are treated 
explicitly. Additional levels up to n = 14 are computed in LTE 
relative to the ground state of Cv. Oscillator strengths from 
ab-initio calculations using the Breit-Pauli /?-matrix method 
(Nahar l2002b) are adopted. Photoionization cross-sections are 
also taken from Nahar's webpage. Effective collision strengths 
for electron impact excitation of transitions among the lowest 24 
fine-structure levels are taken from Aggarwal & Keenan (I2004I I 
and subsequently co-added. All remaining transitions, as well 
as collisional ionization, are treated in analogy to C ii. 

The resulting C ii/m/iv model atom accounts for more than 
1300 radiative and more than 5300 collisional transitions, 
~200 LS' -coupled energy levels and over 20000 frequency 
points - the latter allow the detailed resonance structure of the 
photoionization cross-sections to be sampled with <0.5 A reso- 
lution. Accuracies of the atomic data can range from a typical 
10-20% for ab-initio computations to orders of magnitude for 
approximation formulae. Finally, Voigt profiles are adopted in 
the formal solution using Surface. Wavelengths and oscillator 
strengths of most of the observed transitions are taken from 
Wiese et al. (1996). For Cii ^^6151.3/5 a nd 64 61.9 A the 
wavelengths are adopted from Kurucz & Bell (1 19951 1 and the /- 
values from N02a. Radiative damping parameters are calculated 
from OP lifetimes and coefficients for collisional broadening 
by electron impact are adopted from Griem (119741 for the C ii 
/14267A doublet) or computed according to Cowley (|1971| l. 
Detailed tabulations from quantummechanical computations for 
Stark broadening of several C iv transitions (Schoning I 19931 1 are 
also used. The spectral lines used for abundance analysis are 
listed later in Table 3. 

3. Simultaneous Cii/iii/iv model-atom calibration and 
atmospheric parameter determination 

A reliable set of atomic data considered in model atoms for 
non-LTE calculations can be selected by demanding that the 
model reproduce the observations over the entire stellar parame- 
ter space of relevance. This procedure can be regarded as a 'cali- 
bration' and requires precise atmospheric parameters for the test 
sample stars. These should be free of systematic errors in order 
to prevent one from being mislead when optimising the selection 
of input atomic data. Unfortunately, the parameters of stars are 
a priori not known, they also need to be infeiTed from observa- 
tions. Consequently, a simultaneous solution for an optimal set 
of input atomic data and all stellar parameters is required. These 
parameters should also allow the hydrogen and helium spectra 
to be reproduced. 



^ http://physics.nist.gov/PliysRefData/ASD/indices.litml 
^ http://www-astronomy.mps.ohio-state.edu/~nahar/px.html 



A sample of six apparently slow-rotating B-type dwarfs and 
giants (spectral class B0-B2, luminosity class V-III) randomly 
distributed in Galactic OB associations and in the field of the so- 
lar vicinity are taken as calibration stars. The observational data 
consists of high-resolution spectra with broad wavelength cov- 
erage at very high-S/N (up to ~800 in the Johnson Z?-band), ob- 
tained with FEROS on the ESO 2.2m telescope (La Silla, Chile), 
see Paper 1 for details. 

3.1. Extensive iteration on stellar and atomic variables 

The atmospheric parameter derivation and selection of input 
atomic data are simultaneously performed in an extensive iter- 
ation process. When possible, we try to separate the effects of 
the fundamental parameters and the atomic data on the synthetic 
spectra - the basis for the comparison with observation - in or- 
der to achieve a better understanding of the problem. This is fa- 
cilitated by boundary conditions, like the ionization balance (all 
ionization stages of an element are required to indicate the same 
abundance) or the rules and regularities of atomic physics. 

The iteration is performed on effective temperature Tf^g and 
surface gravity log g, as well as micro-, macroturbulent and pro- 
jected rotational velocities (^, ^ and v sin /, respectively), helium 
and carbon abundances (hereafter e(Hej3 and e(C), respectively) 
and different sets of atomic data. Only the metallicity is fixed 
to a standard solar value (Grevesse & Sauval 1998), a not too 
critical assumption which is furthermore validated a posteriori 
(Przybilla et al. 2007, in prep.). This comprises an enormous 
number of variables (atmospheric parameters and atomic data) 
in the iterative scheme summarised as a flux diagram in Fig.[T] 
The first step concerning the H/He spectrum is solved in Paper 1 . 
There, the He i/ii ionization equilibrium is the main indicator for 
Teff (for the hotter stars), all Balmer lines for log g, the He ii lines 
for ^ and all He lines for ( and v sin /. The second step, involving 
carbon, is required for a fine tuning of the atmospheric parameter 
determination since the metal lines are more sensitive to param- 
eter variations than the H and He lines. Therefore it is possible 
to derive them with a better precision than only from hydrogen 
and helium but at the same time consistently within the error 
limits. Effective temperature and log^ are refined by establish- 
ing the C ii/iii/iv ionization equilibrium in the hottest stars and 
the C ii/iii ionization balance in the cooler stars. The microturbu- 
lent velocity is inferred in the standard way by demanding the 
carbon abundances of the individual lines to be independent of 
equivalent width. Macroturbulent and projected rotational veloc- 
ities are determined by detailed fitting of the carbon line profiles. 
Line fits are performed on the basis of small grids of synthetic 
spectra with different ^ and e(C) via^^-minimisation. 

Since different sets of atomic data are available we have, first 
of all, to intercompare them and judge their reliability to mini- 
mise the uncertainties in the e(C)-determination. Starting from 
an initial model atom we use our expertise for stepwise improve- 
ment, guided by a careful analysis of the atomic structure of the 
model ions and the reactions of the spectrum synthesis to pa- 
rameter variations. Over 100 models were built until consistency 
with all sample stars was achieved simultaneously. Some exam- 
ples of the effects from different model atom realisations on the 
spectrum synthesis are given in the next section. 

The atmospheric parameters derived from C ionization equi- 
librium in Step 2 are verified by re-iterating Step 1 as a final 
check for consistency. Excellent agreement with the observed H 
and He spectra is obtained in all cases, simultaneously in the vi- 

* using the standard logarithmic scale £(X) = log(X/H) + 12 
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Fig. 1. Flux diagram of our extensive iterative procedure for the carbon model atom calibration and the simultaneous precise de- 
termination of stellar parameters, as applied to each programme star, rbb/cbb: radiative/coUisional bound-bound, rbf/cbf: radia- 
tive/collisional bound-free transition data, ew: equivalent width. Step 1 is discussed in Paper 1, the sensitivity of C lines to different 
sets of atomic data in Sect. l3.2l and the sensitivity of C to parameter variations in Sect. l3.4l of the present paper. See the text for details. 



sual and the near-IR (where available). In addition, consistency 
of the observed and model spectral energy distributions (Paper 
2) is found for the final set of stellar parameters. 

By application of the procedure to all programme stars it is 
possible to calibrate the C ii-iv model over the entire parameter 
range (21 500 < T^g < 32 000 K, 3.10 < log^ < 4.30, for dwarfs 
and giants), resulting in a final reference set of atomic data. 
Note that the high-quality spectra available to us are essential 
for this success. They allow us to analyse a wide variety of 
C lines, many of which have never been considered before 
in the study of early B stars. Some weak lines turned out to 
be highly sensitive to non-LTE effects and/or to atmospheric 
parameter variations, namely the Civ lines and CuAA615l and 
6462 A, which change from absorption to emission at higher 
temperatures (discussed in Paper 2). The reproduction of the 
observed trends, despite this high sensitivity, puts strong con- 
straints on the robustness of the final model atom. Note also that 
the presence of lines from three C ionization stages, in hotter 
stars of our sample, allows Teff and \ogg to be derived from the 
ionization equilibrium alone, independent of other indicators. 
The strength of this calibration lies in the simultaneous analysis 
of the large number of C lines of different ionization stages 
in stars covering a wide parameter range. In this way we are 
able to constrain a final set of atomic data independent of any 
specific stellar atmosphere environment. This reference model 
can be used in further applications with a simplified iterative 
scheme, where the only remaining variables are the atmospheric 



parameters and the C abundance. 

3.2. Sensitivity of carbon lines to atomic data 

Reliable level populations are a prerequisite for an accurate non- 
LTE analysis. They can be obtained only when: i) the local tem- 
peratures and particle densities are known (i.e. the atmospheric 
structure); ii) the radiation field is realistic; m) all relevant pro- 
cesses in the statistical equilibrium equations are taken into ac- 
count and iv) high-quality atomic data are available. In partic- 
ular, i) and ii) require a realistic physical model of the stellar 
atmosphere (see Paper 1 for a discussion of the hybrid non-LTE 
approach) and an accurate atmospheric parameter determination. 
Points m) and iv) are related to the model atoms for the non-LTE 
calculations. Shortcomings in any of i)-iv) result in increased 
uncertainties/errors of the analysis. 

Due to the interdependency of all transitions (over 6000 in 
this case) and the non-local character of the radiation field, even 
a restricted non-LTE problem like the one investigated here is 
highly complex. It is impossible to quantify a priori the sensitiv- 
ity of the spectral lines to variations of some of the atomic input 
data. Therefore, one of the few remaining reasons for the large 
spread of carbon abundances found in literature (see Sect.|5]l may 
be different realisations of model atoms (levels/transitions con- 
sidered, atomic data, approximations). Choosing an optimum set 
of input atomic data is not trivial and the construction of reli- 
able model atoms for non-LTE calculations requires a calibra- 
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Fig. 2. Grotrian diagram for the Cii doublet and quartet spin 
systems. The observed multiplets in our spectra analysed here 
are identified. Levels: those marked in bold correspond to levels 
discussed in Figs. |3] and |6l Multiplet transitions: those marked 
by thick lines are highly sensitive to variations of input data for 
photoionization and coUisional cross-sections and for collisional 
excitation (discussed in Figs. |4l|6] and |7] respectively). The lat- 
ter and those marked by thin lines are considered in our linelist 
(Table 3) for abundance derivation. Those marked by dotted lines 
are excluded from the analysis because of contamination with 
telluric lines {AA6251-59, 7231-37 A) or they are too weak even 
at low Tetf (/15889-91 A). Nevertheless they are accounted for in 
the calculations of the level populations. 



tion, guided by extensive comparisons with observation. In the 
following we summarise the experiences made in this process. 

The comparison of observation with model spectra for 
C m and C iv, as computed with our initial model ions, reveals 
little need for improvement in the studied stellar parameter 
range. Both ions are relatively simple, showing (earth)alkali 
electron configurations, which pose little challenge to ab-initio 
computations. However, the actual choice of radiative and col- 
lisional data turned out to be a critical factor for line-formation 
computations of non-LTE-sensitive transitions in C ii, which is 
known to be problematic from the literature (e.g. Lambert [l993l 
Sigut 1 19961 1. This was shortly summarised in Paper 2 but not 
discussed in detail. Here we provide a few comparisons of 
atomic data available from the literature. The influence of 
atomic data on selected Cii lines is also addressed. These 
transitions are highlighted in the Grotrian diagram of our C ii 
model (Fig. |2]l, which will help us to illustrate some of the 
channels leading to the marked non-LTE sensitivity. 

Photoionization cross-sections. The strength of a spectral 
line can be strongly influenced by photoionizations which may 
impact level populations decisively. On the other hand, pho- 
toionization rates depend implicitly on the level populations. 
The largest contribution to the integral Rij comes from frequen- 
cies with large flux and large cross-section. The flux maximum 
in early B stars is located longward of the Lyman jump. The 
ionization of C n is essentially determined by the rates from the 
highly-populated ground state and the low-excitation levels. The 
relative importance of the low-excitation levels may be strength- 
ened in cases where the ground state ionization potential coin- 



cides with that of a major opacity contributor. He i in the present 
case. Photoionizations from the ground state may then be less 
efficient because of the reduced stellar flux shortward of the ion- 
ization edge. Recombinations, on the other hand, are important 
for the population of high-excitation levels (preferentially at high 
i.e. states with large statistical weight), which couple to the 
low-lying states via recombination cascades. For the case of C ii 
this implies: a) an increased sensitivity of the C ii/iii ionization 
balance to the exact run of the photoionization cross-sections of 
levels at low excitation energies and b) an increased sensitivity 
of transitions like Cii .l.i6151, 6462 and 4267 A to non-LTE ef- 
fects because of their participation in the recombination cascade. 

A comparison of total photoionization cross-sections from 
OP (Yan & Seaton [19871 ) and NP97 for two levels (marked 
in the Grotrian diagram. Fig. |2]) is given in Fig. [3] The first 
term, 2p^ ^S, is not directly involved in the formation of the C ii 
/I/} 4267, 6151 and 6462 A transitions. However, it is populated 
considerably and therefore contributes to the Cii/iii ionization 
balance. The cross-section for photoionization from this level to 
the ground state of C iii is small because of radiative selection 
rules. Shortward of ~700A the cross-section rises by about 
two orders of magnitude as the optically allowed channel of 
photoionization to the first excited level of Cm opens. The 
overestimation of the photoionization rate from the 2p^ state 
by adopting the total instead of partial cross-sections is however 
insignificant. The contribution of wavelengths below ~700 A to 
Rij is e.g. less than 1% in the example discussed below. Note the 
wavelength shifts in the resonance structures of the OP and the 
NP97 data. This results in different contributions of the region 
between the threshold and the Lyman edge to the integral in Rij, 
which affect the photoionization rates considerably. On the other 
hand, the photoionization cross-sections for 4f ^F" agree well, 
except for the resonance structure at shortest wavelengths where 
the stellar flux becomes negligible. Consequently, both data (and 
many others for highly-excited levels) are exchangeable without 
showing consequences for the spectrum synthesis computations. 
Note that we adopt experimental threshold wavelengths. On 
the other hand, for 2p^ and other levels at lower excitation 
energy there is a non-negligible effect on the AA4261 and 
6151 A transitions when exchanging both data, as can be seen 
in Fig. |4] In an extreme case, accounting for photoionization 
cross-sections from NP97 for all levels results in a very strong 
Cii AA4261 a line. A reduction of the C abundance by up 
to ~0.8dex is required to fit the observed line profile in our 
calibration stars with such a model atom. As a consequence of 
this we give preference to the OP data over the cross-sections 
of NP97 in our final model atom ('model of reference'), which 
helps us to reproduce observation over the entire parameter 
range consistently. 

Oscillator strengths. Comparisons between multiplet /- 
values from three sources of ab-initio computations are shown 
in Fig. |5] data based on i) the multiconfiguration Hartree-Fock 
method in the Breit-Pauli approximation of FFT04, ii) the Breit- 
Pauli /?-matrix method of N02a and iii) results from /?-matrix 
calculations assuming LS-coupling, as obtained by the Opacity 
Project (Yan et al. I19871 I. The primary source of /-values is 
FFT04, which should be most accurate. Our preference for the 
OP over the N02a data is motivated by the good agreement of 
the former with FFT04 (with two exceptions), while oscillator 
strengths from N02a may show large differences for several 
lines. Data from N02a is therefore adopted only in the cases 
where the other sources do not provide information. 
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Fig. 4. Sensitivity of line profiles of two C ii transitions to vari- 
ations of photoionization cross-sections. The calculations are 
made for three different model atoms using the same model 
atmosphere (as appropriate for HR5285): our reference model 
atom, the reference model atom with cross-sections for 2p^ 
from NP97, and the reference model atom with all photoioniza- 
tion data replaced by values from NP97. The profiles are not 
convolved for effects of rotation or instrumental broadening. 



Collisional ionization cross-sections. Details of the coUi- 
sional ionization cross-sections are important mostly for 
high-excitation levels («>5). Only for these a significant amount 
of electrons in the Maxwell distribution are energetic enough 
to overcome the threshold for the reaction. We use the approxi- 
mation of Seaton i 



() 19621 1. which is expected to be accurate to an 
order of magnitude at best, to evaluate the collisional ionization 
rates because of a lack of any data from ab-initio computations. 
Good agreement between the modelling and most observed 
Cn lines is achieved (see below) when using the standard 
recipe of applying threshold photoionization cross-sections, 
except for two observed transitions. Improvements for the C ii 
/1/16151 and 6462 A lines may be obtained by an empirical 
scaling of threshold cross-sections for the 6f ^F" and 6g 
levels, the upper levels of the transitions. An increase by a 
factor 10, i.e. within the expected uncertainties, helps to repro- 
duce observation over the entire atmospheric parameter range. 
The effects of scaled collision rates for those line profiles are 



-2 
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• ' FFTO-I 



-6 



-2 



Fig. 5. Comparison of oscillator strengths from Nahar (I2002al 
N02a) and Froese Fischer & Tachiev (2004 FFT04) vs. values 
from the Opacity Project (Yan et al. [19871 OP). 



shown in Fig.|6] as well as Cn A4261 A, which remains practi- 
cally unaffected. Larger scaling factors are empirically excluded. 

Collisional excitation cross-sections. The collisional data 
used in a model atom will have an influence on the predicted line 
profiles. Accurate data from ab-initio calculations for larger sets 
of transitions have become available only recently. We employed 
effective collision strengths of Wilson et al. (120051) to construct 
the Cn model. A later revision of part of the data (Wilson et 
al. 120071 ) had negligible influence on the predicted line profiles 
of almost all observable transitions, except for Cii /l/l 39 18/20 
and 6578/82 A, see Fig. |7] for an example. The good agreement 
of abundances derived from these four with other transitions was 
broken when using the improved Wilson et al. (2007) data, re- 
quiring abundance adjustments of up to ~0.3 dex to match ob- 
servation. We find that the situation may be improved for the 
stars of our calibration sample by increasing the effective colli- 
sion strength for the 3s ^S-3p ^P° transition by an empirical fac- 
tor of two, see Fig.|7] This is larger than the typical uncertainty 
of such ab-initio data, which amounts to an estimated 10-20%. 
However, a closer inspection of the energy-dependent collision 
strength for this transition shows that resonances dominate Q. 
in the region near threshold. The positions and strengths of the 
resonances are sensitive to the details of the atomic data calcu- 
lations, in particular to the assumptions made for constructing 
the target. Consequently, more comprehensive ab-initio compu- 



8 



Nieva & Przybilla: Non-LTE line formation for C n-iv in early B-type stars 



C II 4267 




A (A) 

Fig. 6. Reactions of three highly non-LTE-sensitive Hnes to 
changes of colHsional ionization cross-sections. The modifica- 
tions are made for two energy levels which are directly involved 
in the formation of C ii AA6 151 and 6462 A and indirectly in C n 
/14267A. The calculations are made for different model atoms 
with specific values of the reaction cross-section at threshold 
and the same set of atmospheric parameters (as appropriate for 
T Sco). 
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Fig. 7. Effect of employing different effective collision strengths 
on the C ii AA 6578 and 3920 A lines (for HR 1861). Only varia- 
tions of the data from ab-initio calculations (Wilson et al. 120051 
120071 1 are considered. 
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Fig. 8. Comparison of synthetic Cii/iii line profiles for two of 
our sample stars using different model atoms. One calculation 
uses our reference model atom and the other adopts standard 
approximation formulae for evaluating collisional rates for all 
transitions. 
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Fig. 9. Abundance differences from the analysis of CiiAA4267 
and 6578 A (highly sensitive to the input atomic data) and 
Cn/I5145 A ('in LTE'). Displayed are results for our six sam- 
ple stars, as a function of effective temperature. A comparison 
of our initial model atom (all radiative transition data for C ii 
from N02a, NP97) and our final model atom after the calibration 
is shown. See the text for details. 



tations are required to investigate this in detail. However, these 
are beyond the scope of the present work. 

Collisional data from ab-initio calculations are typically 
available only for transitions between relatively low-lying en- 
ergy levels. In the case of C ii the dataset is complete for levels 
up to principal quantum number n - A. Therefore, for the bulk 
of the transitions approximation formulae have to be applied. 
However, trends and regularities from the ab-initio data may be 
used to improve on the standard approximations made for these. 



e.g. dropping the assumption of an energy-independent Q = 1 
for the evaluation of the Allen (119731) formula (see Sect. l2.21 l. In 
practice, simple approximations are used to describe collisional 
processes in most of the model atoms available for non-LTE cal- 
culations at present. One can ask what the effects on synthetic 
profiles will be from such a simplification. This is shown in 
Fig. [8]for selected lines in two of our programme stars. Notable 
corrections in abundance would be required to reproduce the re- 
sults from our reference model, an increase in some of the lines 
but also a decrease in others. Our overall good agreement with 
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observation (Sect.|4| would be destroyed: the non-LTE sensitive 
lines and the 'lines in LTE', which are unaffected by such modi- 
fications, would indicate widely different abundances. Note that 
the effects vary from star to star. 

We conclude the discussion on the impact of atomic data 
on line-formation calculations in Fig. |9] Here, we quantify 
the differences in abundance derived with our initial model 
atom (built from the available homogeneous set of atomic 
data, NP97/N02a) and with our model of reference (after the 
calibration) for the entire sample of stars. Abundances from 
the strong and non-LTE-sensitive Cn AA4261 and 6578 A 
transitions are compared to abundances derived from the weaker 
Cn /15145A line, which is almost insensitive to the details of 
the model atom (i.e. it is 'in LTE'). The relative abundance is 
displayed as a function of effective temperature for each star 
(see Section ID for final results). Pronounced systematic trends 
exist when the initial model atom is used, correlating with the 
strength of the lines (i.e. stronger lines show a larger sensitivity 
to non-LTE effects and therefore to the input atomic data). 
These trends and abundance differences from C n AA4261 and 
6578 A almost vanish when our model of reference is used. The 
remaining small differences may be reduced even further when 
improved atomic data, in particular for collisions involving 
high-excitation levels, become available. Note, that the Cii 
/16582 A line (not displayed here) shows a trend similar to that 
of /16578A. The whole multiplet centred on Cii /15145A also 
behaves consistently. This kind of test has been made for all 
C ii-iv lines for every set of input atomic data in the empirical 
calibration process of the model atom. The detailed study of 
non-LTE and LTE abundances in Sect. |4] helps to identify the 
lines in LTE for the different atmospheric parameters under 
analysis. These may provide good starting points for further 
analyses when one desires to avoid non-LTE effects. 

3.3. Line-formation details 

A closer study of the underlying line-formation processes allows 
the nature of the non-LTE effects to be understood. This is shown 
in Fig.[TO]for some representative transitions of Cii/iii/iv in our 
programme star r Sco. Departure coefficients = nj^^™/nj"™ 
are displayed in the left panel of Fig. [10] for the levels involved 
in the transitions of interest and the ion ground states. Non-LTE 
departures of the level occupations impact the line source func- 
tion S\. We recall that the ratio of 5 1 to the Planck function By 
(shown for selected transitions in the right panel of Fig. [TOl l is 
Si/By = [expihvij/kT) - l]/[bi/bjexpihvij/kT) - 1]. Si/By is 
determined by the ratio of the departure coefficients for the lower 
and upper levels (/, j). An overpopulation of the upper level rel- 
ative to the lower (i.e. Si/By > 1) results in non-LTE weaken- 
ing of the line and may lead to emission in cases of pronounced 
overpopulation, while the inverse gives non-LTE strengthening. 

The level populations reach detailed equilibrium values deep 
in the atmosphere, where large collisional rates and small mean- 
free paths between photon absorptions (both because of the high 
densities of the plasma) enforce this inner boundary condition. 
Double-ionized carbon is the main ionization stage at the tem- 
peratures of T Sco - the C III ground state is close to LTE. Single- 
ionized carbon is overionized at line-formation depths, therefore 
the levels are underpopulated relative to LTE, and C iv and the 
C V ground state are overpopulated. 

In general, level populations in C ii depart most from detailed 
equilibrium in the low-excitation states and approach LTE values 
gradually with increasing excitation energy, as collisions facili- 
tate coupling with the C iii ground state. Therefore, most of the 



non-LTE-sensitive transitions in C ii have upper levels that are 
overpopulated relative to the lower level, such that the lines ex- 
perience slight (A 3920 A) to notable weakening (A 4267 A) rela- 
tive to LTE and may even turn into emission (AA 6151, 6462 A). 
The /l/l 6578/82 A doublet experiences non-LTE strengthening 
for lower effective temperatures, while in hotter objects like 
T Sco the lines are found to be close to LTE. The other observ- 
able C II lines arise in the quartet spin system. They are weaker 
(C II /1 5 145 A is the strongest), i.e. they are formed deeper in the 
atmosphere, and their formation involves high-excitation levels 
which are coupled collisionally at these depths, such that these 
lines are essentially in LTE. Note that the behaviour of Si/By for 
A 4267 A is in good agreement with the findings of Sigut ( 1 19961 1. 
On the other hand, notable differences exist for the AA 6578/82 A 
doublet in particular at higher T^ff. The reasons for this will be 
discussed in the next section (but see also Fig.[T4b. 

The Cm transitions can also experience both, non-LTE 
weakening (like the strong /J4187A) and non-LTE strength- 
ening (like the strong triplet 4647-4651 A). The Civ doublet 
AA 5801/12, which becomes observable only in our hottest stars, 
shows a pronounced non-LTE strengthening. 

3.4. Sensitivity of e(C) to atmosplieric parameter variations 

Spectral lines of carbon, like many other metal lines, can react 
sensitively to variations of the stellar atmospheric parameters. 
This property provides us with a powerful tool for the atmo- 
spheric parameter and abundance determination, using our itera- 
tion scheme. This is a non-trivial and crucial step in the analysis, 
which must be performed carefully in order to avoid system- 
atic error. We investigate consequences of systematically biased 
atmospheric parameters on carbon line profiles and the derived 
abundances in the following 0. 

The offsets for the parameters (in effective temperature 
A Tsff- -2000 K, surface gravity A log g= +0.2 dex and micro- 
turbulent velocity A ^= +5 kms^') are representative for system- 
atic discrepancies between our final values and those from pre- 
vious studies (and also among previous studies). Note that they 
are much larger than our statistical uncertainties. Such discrep- 
ancies may be caused by several factors, among others: i) photo- 
metric effective temperature calibrations based on model atmo- 
spheres with insufficient line blanketing, ii) spectroscopic ioniza- 
tion equilibria based on predictions of incomplete model atoms, 
III) the assumption of LTE for the computation of Balmer line 
profiles, which may cause a systematic overestimate of log g by 
up to 0.2 dex (see Paper 1). 

Figure [TT] shows the sensitivity of selected Cii/iii/iv lines in 
two giants of our sample to variations of like Teff , log g and ^. 
The profiles accounting for these variations are compared with 
those computed with our final atmospheric parameters and an 
averaged carbon abundance of e(C) = 8.32. The sensitivity to pa- 
rameter variations differs from line to line. 

For the hotter stars, the C iv and - when strong enough - 
the Cii/16151A multiplets are ideal indicators for both Tetf and 
log g, while the rest of the lines are mostly sensitive to changes in 

' Note that an underestimated source of systematic error in view of 
line blanketing are abundance 'standards' in the the widely used clas- 
sical Kurucz (1993a) ODFs, computed on the basis of (scaled) solar 
abundances according to Anders & Grevesse ( 1989). Employing ODFs 
with appropriately reduced metallicity for the model atmosphere calcu- 
lations (see Paper 1 for a discussion) decreases the atmospheric tem- 
perature structure via a backwarming effect, by up to ~500K in the 
line-formation region, i.e. slightly higher than our uncertainties in T^fi. 
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Fig. 11. Sensitivity of selected C ii/iii/iv lines to atmospheric pa- 
rameter variations in two giants: HR3055 (hotter, BO III) and 
HR3468 (cooler, B 1,5 III), Our solution corresponds to the fi- 
nal parameters from Table |2] and a constant value of e(C) for all 
lines. The parameter offsets are typical for statistical and system- 
atic uncertainties from published values. The theoretical spectra 
are convolved with a rotational profile for vsini - 20 km s 
Our solution establishes the ionization equilibrium. The star 
HR 3468 is too cool to show C iv lines. 



Table 1. Systematic uncertainties in carbon abundances (in dex, 
relative to our final results. Table 3) caused by atmospheric 
parameter variations and the assumption of LTE for the line- 
formation calculations. 
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Tef[ and the strong C ii/m lines also to changes in ^, For the cooler 
stars, the Cii AA 4267 and 6151 A multiplets are highly sensitive 
to variations of T^-ff and ^, the rest of the lines are mostly sensi- 
tive to Teff and the strong lines also to ^, 

The response to variations in log g can be amplified for C ii 
AA 6578/82 A at higher gravities, because they are formed on the 
red wing of He The local continuum and therefore the line- 
formation depths may change as a consequence, impacting the 
non-LTE effects, A correct treatment of the hydrogen B aimer 
line-opacity therefore plays an important role in this context. 
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Fig. 12. Sensitivity of carbon abundances 
to stellar parameter variations: T^ff (up- 
per panel), log g (mid panel) and mi- 
croturbulent velocity (lower panel). The 
offsets of the parameters are displayed 
in the lower left part of each panel. A 
large spread in abundance from individ- 
ual lines of the three ionization stages re- 
sults, in particular for variations of Teff 
(upper panel). Note also the implications 
of using only few lines of one ionization 
stage for abundance determinations in the 
presence of systematic errors in the atmo- 
spheric parameters. The grey bands cor- 
respond to Icr-uncertainties of the stellar 
carbon abundance in our solution (sum- 
marised in Table |2j. 



Non-LTE effects strengthen the Balmer line wings in particular 
at higher temperatures (Paper 1). 

For HR 3055 we have quantified the systematic effects ex- 
emplified in Fig. [TT] (note that C iii AA6A1 A is blended and C ii 
/16151 A is too weak in this case) by deriving non-LTE carbon 
abundances for the modified values of Tefj , log g and ^ and com- 
paring them to our solution for individual lines. This is sum- 
marised in Table[T]and visualised in Fig. [12] Table[T]also shows 
systematic offsets that arise from the assumption of LTE for the 
line-formation calculations. This comparison allows us to iden- 
tify the relative importance of atmospheric parameters/non-LTE 
effects for some key spectral lines. Note that the C iv lines are ex- 
tremely sensitive to changes in T^ff and log g at this temperature 
(31200 K), with discrepancies amounting to up to ~+1.0dex 
in abundance for Areff= -2000 kQ. The Cii/m ionization bal- 
ance is also never established for a variation of these parame- 
ters (abundance changes down to -0.40 dex for C ii and up to 
-1-0.35 dex for C m when compared to our solution). An expected 
reduction of the abundances from strong lines is obtained for an 
increased microturbulence. Note that the systematic variations of 
carbon abundance with ^ for some lines are significant consid- 
ering the high accuracy we are aiming at, despite smaller effects 
in general than for Teff and log g variations. 

The solutions for the modified atmospheric parameters are 
characterised by a large scatter of C abundances from indivi- 
dual lines (the statistical 1-cr uncertainties increase by up to ~0.4 
dex). Note that variations of Teff show by far the largest effects. 

We conclude that the use of only a few spectral lines from 
one ionization stage for C abundance determinations - which is 
common practice in the literature - may have serious implica- 
tions on the final accuracy. Possible systematic discrepancies, as 
indicated here, may remain unrecognised. Moreover, the oppor- 
tunity to improve on the atmospheric parameter determination 
by establishing the highly parameter-sensitive ionization balance 
may be missed in such cases. This high sensitivity of the carbon 
lines to atmospheric parameter variations can be used as an im- 
portant tool for precise quantitative spectral analyses of this kind 
of star when applying the calibrated model atom in the future. 

Some lines are practically unaffected by non-LTE effects, as 
indicated by Table [1] (see also Sect. 13.2b . This implies that they 



^ Even larger offsets in T^g are found with respect to the literature, 
up to AT^ff ^ -4000 K, see Section 5. 



are almost insensitive to any reasonable choice of model atom. 
On the other hand, they are highly sensitive to the choice of 
atmospheric parameters. This property is highly useful for the 
model atom calibration because it helps to disentangle effects 
due to non-LTE from those due to inaccuracies in the stellar pa- 
rameters, facilitating a reduction of systematic errors. Typical 
systematic uncertainties in the atmospheric parameters can have 
a similar or even larger impact on the carbon abundance deter- 
mination than a neglect of non-LTE effects. 

4. Results 

Accurate atmospheric parameters and carbon abundances are de- 
rived from line-profile fitting by minimisation, which puts 
tighter constraints than matching only equivalent widths. These 
are summarised in Table |2] for our six programme stars (as ob- 
tained from our iterative process. Fig. [1]). The uncertainties of 
Teff and \ogg are estimated from the extremely sensitive car- 
bon ionization equilibrium. They are lower for hotter stars be- 
cause of the additional restrictions imposed by the presence of 
lines from three ionization stages. For velocities, we provide the 
Icr-uncertainties from the analysis of the entire line ensemble. 
Projected rotational velocity and (radial-tangential) macroturbu- 
lence (Gray [19921 p. 407ff.) were simultaneously derived allow- 
ing for small line-to-line variations in order to obtain an opti- 
mum fit (see Ryans et al. 120021 for the case of B supergiants). 

We should emphasise that the final set of atmospheric pa- 
rameters as derived from the C ii/iii/iv ionization balance is in 
agreement with those from our previous quantitative analyses of 
these stars. In particular, a simultaneous match is achieved for 

i) the H and He lines in the visual and (where available) in the 
near-IR, including the He i/ii ionization equilibrium in the hot- 
ter stars (Paper 1, also the source of helium abundance y) and 

ii) the spectral energy distributions from the UV to the near IR 
(Paper 2). A similar degree of consistency is typically not ob- 
tained in comparable studies of early-type stars. 

A large quantity of carbon lines is analysed for the first time, 
giving consistent abundances for all of them. The excellent qual- 
ity of the observed spectra combined with our improved analy- 
sis technique allow us to achieve such precise results. The Icr- 
uncertainties from the line-to-line scatter are typically of the or- 
der 0.05-0. 10 dex. We estimate the systematic uncertainties due 
to remaining eiTors in atmospheric parameters and atomic data 
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Fig. 13. Non-LTE and LTE abundances de- 
rived from line profile fits to individual Cn- 
IV lines in the programme stars as a function 
of equivalent width. The ID of the stars and 
the spectral classification is given in the up- 
per left comer of each panel. The grey rectan- 
gles correspond to Icr-uncertainties of the stel- 
lar carbon abundance from the line-to-line scat- 
ter. Identification of lines with high sensitiv- 
ity to non-LTE effects is displayed. Emission 
lines are marked by crosses (Cii /1/16151 and 
6462 A in T Sco and /16462 A in HR 1 86 1 ) : LTE 
calculations are not able to reproduce them 
even qualitatively. 



Table 2. Stellar parameters and carbon abundances of the programme stars. 



tSco 
HD 149438 



HR3055 
HD 63922 



HR1861 
HD 36591 



HR2928 
HD 61068 



HR5285 
HD 122980 



HR3468 
HD 74575 



ScoCen Field OriOBl Field ScoCen Field 

32000 ±300 31200 ±300 27000 ± 400 26300 ± 400 21500 ±400 22900 ±400 



Teff (K) 
logg(cgs) 
f (km s-') 
vsini(km s"') 
^(km s-') 



4.12 i 

34 

12: 



4.30 ± 0.05 3.95 ± 0.05 

5± 1 8± 1 

4 ± 2 29 ± 4 

4 ± 2 37 ± 8 

y (by number) 0.089 ± 0.01 0.080 ±0.01 0.089 

£(C)(dex) 8.30 ±0.12 8.33 ± 0.08 8.33 d 

#C lines 32 19 30 



0.05 
1 

: 1 



4.15 i 

3d 

14: 
20: 

:0.01 0.089: 

0.08 8.27 4 



0.05 
1 

: 1 
:2 

:0.01 

0.07 



22 



4.20 ± 0.05 
4± 1 
18 ± 1 

0.089 ±0.01 
8.29 ± 0.05 
22 



3.60 ± 0.05 
5 ± 1 
11 ±2 
20 ± 1 
0.089 ± 0.01 
8.37 ±0.10 
19 



to be of the order 0.10-0. 15 dex, taking Table [T] and our experi- 
ences from Sect. 13. 21 as a guideline. 

Details on the analysis of individual lines can be found in 
Table 3. This summarises line identifications, level designations, 
log gf values, excitation potential of the lower level xi, and for 
each star equivalent widths Wa and the derived non-LTE and 
LTE abundances. Note that we analyse as many lines as possi- 
ble per star, excluding only features with strong blends by other 
chemical species. This helps us to understand the behaviour of 
each line and the quality of our modelling at different temper- 
atures and gravities (see the previous section). Transitions in- 
volving autoionizing states are not accounted for explicitly in the 
spectrum synthesis (however, they are considered via resonances 
in the photoionization cross-sections). These are lines like C 11 
AA4075 (a blend with On), 4318 (a blend with Sii), 4374.3, 
4375.1, 4376.6, 4411.2/5 and 4627.4 A, which have sometimes 
been used for abundance determinations in previous studies. 

The high quality of our line fits to the observed spectrum of 
tSco is demonstrated in Fig.[T7] available online, for almost all 
the analysed transitions. Similar information for the other pro- 
gramme stars can be found in Figs[T8]^22]in the online version. 
Of central importance is that the abundances derived from the in- 
dividual lines show a small scatter in each sample star Good fits 
to individual lines can almost always be obtained, however this 
does not imply consistency in the entire analysis. An example 
are the abundances from Fig. [12] which were also derived from 
high-quality line fits, but which show inconsistencies, expressed 



as a large line-to-line scatter, nonetheless. Another example from 
many other studies are discrepant lines that are excluded from 
the analysis in order to reduce the statistical uncertainties. The 
present work improves on this because the underlying physics is 
solved in a more consistent way. 

We find a single exception to our overall good line fits: the 
doublet C II /1/14267.0/2 A in r Sco, where the fine-structure com- 
ponents are resolved (Fig. [TtI i. The synthetic profile is slightly 
broader than the observed one, even when neglecting micro- and 
macroturbulent broadening. Such a detail cannot be observed in 
the other stars because of their higher v sin /. 

The data from Table 3 are visualised in Fig. [13] where abun- 
dances are shown as a function of Wa. Excellent consistency is 
found in the case of the non-LTE analysis while the quality of 
the results is considerably degraded in LTE. Note that the sub- 
set of the strong lines provides non-LTE abundances that are in 
good agreement with those derived from the entire sample of 
lines. This proves that our reference model atom is suited well 
for analyses of low-S/N spectra and fast-rotating stars, where 
only the strongest lines are measurable and where such a consis- 
tency check is not feasible. Note also that the non-LTE effects 
differ from line to line and also from star to star 

As the carbon abundances presented here are derived via 
line-profile fitting of our synthetic spectra, they do not depend 
on equivalent width measurements. However, careful determina- 
tions of Wa via direct integration, cross-checked by Voigt-profile 
fitting, were performed to allow a visualisation like in Fig.[T3] 



Table 3. Carbon abundance analysis for the programme stars. 



tSco HR3055 HR1861 HR2928 HR5285 HR3468 

Ion A (A) l-u log gf Xi (eV) Wx NLTE LTE Wx NLTE LTE Wx NLTE LTE Wx NLTE LTE Wx NLTE LTE Wx NLTE LTE 
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bl: blend with other metal line; ... : too weak or absent; XX: LTE does not reproduce the emission; 

logg/-values in formal solution: Wiese et al. (1996, uncertainties within: 3% A, 10% B, 25% C) and Nahar (2002a, uncertainties within 10-25% N) 
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Comparisons between both methods can give discrepancies 
up to 10-20% in the worst case, mostly for weak lines in spec- 
tral regions where the continuum is difficult to be set, the S/N is 
slightly degraded or there are partial blends with other lines. 

One case requires further discussion: the C ii AA 6578/82 A 
doublet in rSco. The remaining problems may be an artifact 
from the data reduction. The location of the red wing of Ha 
coincides with a bad column of the CCD of FEROS (the spec- 
tral orders are oriented along columns in this spectrograph). A 
perfect correction for this cannot be provided. Different wave- 
lengths are affected, depending on the radial velocity shifts of 
the object. Note that the region around C ii AA 6578/82 A doublet 
has been normalised relative to the local continuum in Fig.fTTl 

Finally, even when non-LTE-insensitive lines are considered 
in the analysis, it is mandatory to verify that the ionization bal- 
ance has been established. This is in order to avoid systematic er- 
rors in the atmospheric parameters (in particular Tgs) and there- 
fore in abundance. It is recommended that such lines in LTE be 
used for abundance analyses when an optimised model atom is 
not available and the S/N of the observed spectra is high enough 
to measure them. Table 3 may act as a guideline to identify ap- 
propriate transitions, depending on atmospheric parameters. 

5. Comparison with previous work 

We concentrate here on the main sources of systematic un- 
certainties that can bias abundance analyses, as identified in 
Sects. 13.21 & 13.41 atomic data, i.e. different model atoms, and 
stellar parameters, i.e. in particular T^ff scales. A more compre- 
hensive comparison is too complex to be performed because of 
the large number of variables involved, which in some cases are 
not even documented. This includes numerous factors related to 
observation (e.g. quality of the analysed spectra, continuum rec- 
tification, equivalent width measurements, line blends), model 
atmospheres (e.g. codes, abundance standards and linelists con- 
sidered for line blanketing) and line-formation calculations (e.g. 
handling of line blocking, oscillator strengths, line broadening). 
A comparison of our results on carbon abundances in the solar 
neighbourhood with previous studies will be made in an extra 
section, because of the wider implications. 

5.1. Predictions from different non-LTE modei atoms 

We have computed non-LTE line profiles based on the present 
model atom of reference and the Eber & Butler C1988. EB88) 
model and have compared them with Sigut's (I1996I I non-LTE 
data and our LTE results in Fig. [14] All models are calcu- 
lated for the same set of atmospheric parameters and carbon 
abundances. Note that only one fundamental difference exists: 
our computations account for non-LTE populations for hydro- 
gen while Sigut assumes LTE. This is one of the reasons why 
Sigut ( 1996) found problems in reproducing observed trends 
for the C ii AA 6578/82 A doublet at T^ff > 25 000 K: Ha was as- 
sumed to be formed in LTE in that work. 

We obtain consistent abundances from the application of our 
model atom of reference to observations for the Cii A 4267 A 
transition, as shown in Sect. H] Consequently, the good agree- 
ment with Sigut's predictions indicates that his model atom is 
also highly useful for abundance determinations from this line. 
Note that Sigut compares his results (for e(C) = 8.55) only qual- 
itatively with observed equivalent widths from different sources, 
using stellar parameters derived in one of that studies, or from a 
photometric Teff-calibration. A reduction of carbon abundance 
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Fig. 14. Predicted equivalent widths of some C ii lines from dif- 
ferent approaches as a function of T^ff. All calculations were per- 
formed for log g - 4.0, ^ =5 km s"' and £(C) = 8.55, in order to 
facilitate a comparison with Sigut ( 119961 . 



and an improved determination of atmospheric parameters as 
proposed in the present work may bring observation in much 
better agreement with his predictions (his Fig. 1). In the region 
~22 000 K < Teff < 28 000 K good agreement is also found with 
the predictions of the EB88 model atom. However, this model 
predicts too large equivalent widths and therefore lower abun- 
dances outside this region. 

Our analysis gives abundances from the /l/i 6578/82 A 
doublet consistent with those from other lines (tSco may be 
an exception, as discussed above). The trends predicted with 
our model atom of reference differ from the three other model 
calculations. Sigut's and EB88's non-LTE predictions agree 
qualitatively up to T^ff ~ 22 000 K, Sigut's model indicating 
larger and EB88 smaller equivalent widths, respectively. 
However, our model shows non-LTE strengthening throughout, 
approaching LTE at the high Teff-limit, while the two other mod- 
els imply pronounced non-LTE weakening at T^^fi > 25 000 K. 
Note that this doublet experiences significant non-LTE effects 
only when the lines are strong, e.g. for high values of carbon 
abundance such as those discussed in Fig. [14] The non-LTE 
effects are abundance-dependent, they are reduced with decreas- 
ing abundance. Finally, the predicted equivalent widths for the 
widely used Cii multiplet including /15145A are practically 
independent of the model assumptions. The lines are close to 
LTE. Any inconsistencies of carbon abundances based on this 
multiplet have to be related to other effects, but not by the choice 
of the model atom. This multiplet, of the quartet spin system, 
was not analysed by Sigut. 
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Fig. 15. Comparison of our final Teft for the sample stars with 
values derived from photometric calibrations: TeffCQ) (Daflon 
et al. [19991 Lyubi mkov e t al. 2002 LRRL) and T,ff([u - b]) 
(Napiwotzki et al. 119931 1. Photometric indices are computed 
from SIMBAD data. Spectroscopic temperatures according to 
Kilian ( 119921 1 are also displayed. 



5.2. Effective temperatures 

Abundance determinations can be systematically biased even 
in when a realistic non-LTE model atom is available. The at- 
mospheric parameters need to be derived accurately as well. 
We have quantified the effects of inappropriately chosen atmo- 
spheric parameters in Sect. 13.41 Largely discrepant abundances 
are found from different ionization stages in particular for inac- 
curate effective temperatures. We were confronted with this from 
the beginning of our work, when adopting effective temperatures 
derived from common photometric calibrations. It was not pos- 
sible to establish the ionization equilibrium of C n/m even for the 
lines in LTE in that case. This motivated us to use an iterative ap- 
proach to constrain all variables that are involved in the analysis 
simultaneously and in a self-consistent way. 

The final values for T^ff are compared with other derivations 
in Fig. [15] The comparison concentrates on widely used photo- 
metric Teff-calibrations, which provide a fast determination of at- 
mospheric parameters, allowing the analysis of larger samples of 
stars to be performed. This includes methods using broad-band 
Johnson photometry, like the reddening-free Q index (Daflon 
et al. [T999i: Lyubimkov et al. |20021 l or small-band Stromgren 
photometry, as the Tgff{[u - b]) calibration of Napiwotzki et 
al. ( 119931 ). The required Johnson and Stromgren magnitudes for 
the individual stars are adopted from the SIMBAD database. 

The discrepancies in Teff can be large, amounting up to 
~4000 K, much larger than the offsets studied in Sect. 13.41 The 
differences in abundance from Cn and Cm can then achieve 
~0.7 dex, similar to what is found in previous studies. Our Teff- 
values are typically higher than those from the photometric esti- 
mates. The present sample is too small to facilitate an improved 
empirical Teff (2)-calibration on its own. However, there is an in- 
dication that giants have cooler atmospheres than dwarf stars at 
the same photometric index. This is qualitatively in agreement 
with the findings of Lyubimkov et al. ( 120021 1. 
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Fig. 16. Comparison of C abundances derived from our sample 
of early B-type stars in the solar vicinity (from the field and from 
OB associations) with results from the literature. Similar objects 
at distances shorter than ~ 1 kpc from the Sun and at galactocen- 
tric distances within up to ~500 pc difference with respect to the 
location of the Sun are considered. The most recent solar values 
are also indicated. The binsize of the histograms is related to the 
statistical uncertainty of each sample. Our sample is not large 
enough for a statistical comparison, therefore a column with a 
thickness of our uncertainty is displayed. The programme stars 
coincide with six objects from Kilian ( 119921 1: her abundances 
show a much larger spread for these. See the text for details 



Spectroscopic Teif-determinations by Kilian (1992) are also 
shown in Fig. [15] The approach is conceptually the most similar 
to the present study despite considerable progress made over the 
past 15 years and, as a consequence, gives the closest agreement 
with our values (except for the cool giant). 

These comparisons show that the effective temperatures de- 
rived from available photometric calibrations in early B-type 
stars should be considered only as estimates. If adopted as final 
temperatures, one should be aware of the potential systematic 
effects on the abundance determination (see Sect. 13.41 ). Instead, 
these could be used as starting points for further refinements 
via well-understood spectroscopic indicators. This is in order to 
bring the model into agreement with all available indicators si- 
multaneously: e.g. multiple H and He lines and metal ionization 
equilibria, as in the present approach. All the atmospheric pa- 
rameters could then be tightly constrained. Such a step is essen- 
tial if the Icr-uncertainties in the abundance analysis are required 
to be smaller than ~0.3 dex. 

6. The stellar present-day C abundance in the solar 
neighbourhood 

Early-type stars can act as tracers for the present-day chemi- 
cal composition in the solar neighbourhood. The present small 
sample of stars is suited to address this topic, as the objects are 
randomly distributed over nearby OB associations (three stars 
in Ori OB 1 and Sco-Cen, see Table [2]i and the field (the other 
stars). They are located at distances shorter than ~1 kpc from the 
Sun and at galactocentric distances within up to ~500pc differ- 
ence with respect to the location of the Sun. We find a highly- 
homogeneous carbon abundance of e(C) = 8.32+0.04 (Icr sta- 
tistical uncertainty). Continuing the discussion from the last sec- 
tion, the carbon abundance derived in this work is compared with 
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Table 4. Carbon abundances of different objects in the solar 
vicinity 



Objects 


b(C) 


Source 


B stars (pristine value) 


8.35+0.05 


present work 


B stars 


8.25+0.08 


Herrero (|2003k 


Orion H ii (gas) 


8.42+0.02 


Esteban et al. t20()4l 


Orion H ii (gas+dust) 


8.52+0.02 


Esteban et al. (2004) 


young F and G stars 


8.55+0.10 


Sotia & Meyer (200n 


ISM 


8.15+0.06 


Sotia & Meyer ( 2001) 


Sun 


8.39+0.05 


Asplund et al. ( 2005 ) 


Sun 


8.52+0.06 


Grevesse & Sauval (,1998J 



results from selected previous studies of early B-type dwarf and 
giant stars in non-LTE, as shown in Fig. [16] 

Non-LTE analyses based on the EB88 model atom, using dif- 
ferent linelists, find systematically higher abundances than pre- 
vious LTE studies. However, a large spread in abundance is still 
present in each sample. Our results can be directly compared 
with those of the purely spectroscopic study of Kilian ( |1992| l. A 
mean e(C) - 8.19+0.12 is derived from the same six stars when 
adopting her results. This comparison is highly important as it 
shows that a self-consistent analysis, like the one presented here, 
may drastically reduce the statistical scatter and may also imply 
a considerable systematic shift of the mean abundance. Similar 
to Kilian, Gummersbach et al. ( 119981 not shown in Fig.[T6]l also 
performed a consistent spectroscopic analysis of five stars in the 
solar vicinity as part of their Galactic abundance gradient study. 
Their results also indicate systematically lower abundances and 
a considerable spread. 

Examples of non-LTE analyses for larger samples of stars 
based on effective temperatures derived from photometric cali- 
brations are also displayed in Fig. [16] Gies & Lambert ( |1992| l 
and Cunha & Lambert ( |1994| l derived Tetf from calibrations 
based on Stromgren photometry and non-LTE abundances from 
the Cii multiplet around /i5145A and where possible from Cii 
/1/15648/62A. Both multiplets originate in the quartet spin sys- 
tem and they are almost unaffected by non-LTE effects. For 
the comparison we have excluded the five supergiant stars of 
the sample of Gies & Lambert ( 1992). Gies & Lambert ( |1992| l 
have two stars in common with the present sample, indicating 
systematically lower abundances: for HR 2928 A£(C) - -0.3 dex 
and for HR 1861 Ae(C) - -0.2 dex (from the same lines). They 
also consider the doublet C ii /1/16578/82 A for several stars of 
their sample, but obtain even lower abundances compared to 
the lines from the quartet spin system (by approximately an- 
other -0.3 dex). Cunha & Lambert (1994) have only one star 
in common with us, HR 1861. They found a similar abundance 
for this star as Gies & Lambert ( 1992). More recently, Daflon et 
al. ( 119991 12001 al l derived temperatures from their calibration of 
the Johnson Q-parameter and non-LTE abundances from the C ii 
multiplet around A5 145 A. 

All the data on carbon abundances from early-type stars sum- 
marised above was derived from interpretation of observation, 
which may be affected by many sources of systematic error (see 
Sect. [3]l. These data are interpreted in turn to test models of 
massive star evolution and the chemical evolution of the Galaxy, 
which are anchor points for studies of stellar and galactochemi- 
cal evolution in general. Any misinterpretation in the first step in 
this chain may have severe consequences for our understanding 
of the whole picture. In this context, it is important to investigate 
the conclusions drawn from a diagram like Fig. [16] 



Previous non-LTE studies all show a broad range of carbon 
abundances in the solar neighbourhood, spanning a factor of ~10 
in total. Current models of stellar evolution and chemical evolu- 
tion of the Galaxy do not provide a physical explanation to the 
large range shown in Fig. [16] for this kind of star located in the 
solar vicinity(3. If this extended distribution were real, its physi- 
cal interpretation may require to account for the following physi- 
cal processes, i) Mechanisms that alter the atmospheric structure, 
such as sufficiently strong magnetic fields or diffusion/radiative 
levitation processes acting on short timescales and independent 
of atmospheric parameters and evolutionary age. These are so far 
discussed only for a few cases, well motivated by observational 
evidence. A common occurrence of these processes would im- 
ply that the basic assumption of homogeneity of the atmosphere 
may be invalid and as a consequence the modelling techniques 
applied so far may be inadequate, ii) Extremely efficient deple- 
tion mechanisms in the course of stellar evolution, by one to two 
orders of magnitude larger than currently predicted (e.g. Meynet 
& Maeder I2003I I. Already on the main sequence mixing with 
matter from the stellar core would require a higher efficiency 
than currently predicted for convection during the first dredge- 
up. Eventually, also enrichment mechanisms may be required 
which contradict nucleosynthesis (the 3Q'-process is not active 
in OB dwarfs and giants, carbon is instead depleted by the CNO 
cycle). Ill) An enormous chemical inhomogeneity of the present- 
day interstellar material in the solar neighbourhood out of which 
the stars have been formed. This would require the possibility to 
change abundances practically instantly even within single clus- 
ters by amounts that are otherwise attributed to the past ~12 Gyrs 
of Galactochemical evolution (e.g. Chiappini et al. 120031) . 

Concerning the chemical inhomogeneity, non-uniform abun- 
dance distributions as a result of galactochemical evolution are 
also discussed in literature. E.g., Oey (2003) indicates a scatter 
of 0. 13 dex in oxygen abundance for a population at solar metal- 
licity, which is in reasonable agreement with previous observa- 
tional constraints. However, this is an upper limit, as homogeni- 
sation processes were neglected entirely, such that the true scat- 
ter can be expected to be much lower. 

The present work, in combination with Paper 1 and 2, avoids 
these fundamental problems. It implies homogeneous abun- 
dances after improving the modelling and the analysis methodol- 
ogy, bringing all model aspects into agreement with observation 
at once, which was not achieved so far. 

i) It conforms with the finding of a uniform gas-phase carbon 
abundance in the ISM (out to distances of 1.5 kpc, e.g. Sofia & 
Meyer 1200 11 and references therein), with the systematic offset 
between the two values being a consequence of dust formation in 
the ISM. The uncertainties of the ISM abundance (see Table ]4]i 
are determined by the accuracy to which the oscillator strengths 
of the resonance lines in the UV are known, and should not ex- 
ceed more than ~10%. 

n) It also agrees with Galactochemical evolution models, 
which predict homogeneous abundances in the solar neighbour- 
hood (e.g. Chiappini et al. 2003). A variety of hydrodynamic 
processes should keep the ISM chemically well-mixed on small 
time-scales (Edmunds 1 19751 Roy & Kunth [T995] l. The variation 
due to the Galactic abundance gradient should amount to up to 
~0.04-0.08dexkpc ', which is of the order of the uncertainty 
found here (one kpc is the maximum Galactocentric distance 
sampled in our comparison). 



' When considering also LTE results (e.g. Kane et al. |1980t the whole 
abundance range rises up to 1.6 dex (a factor of 40). 
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m) Another aspect concerns the notably sub-solar abun- 
dances from early-type stars found so far (e.g. Herrero 120031 
with respect to the old solar standard, Grevesse & Sauval 1998 
- see Table|4]l. The present findings remedy the situation, in par- 
ticular if accounting for the recently revised solar carbon abun- 
dance (Asplund et al. 120051 see also Table|4]l. 

iv) Finally, rotationally-induced mixing with CN-processed 
material from the core may change the atmospheric compo- 
sition. Theory predicts a depletion of carbon by ~0.03dex 
for a star of 20 M© with initial rotational velocity 300km s"' 
evolving from the zero-age Main Sequence to the end of the 
Main Sequence stage, reaching ~0.15dex in the supergiant 
stage (e.g. Meynet & Maeder 2003 ). The presence of a magnetic 
field may amplify rotational mixing (Maeder & Mevnet [20051) . 
however the effects on abundances are not expected to exceed 
a factor ~2. As absolute rotational velocities of stars can be 
measured only for a few exceptionally cases, this topic can be 
addressed comprehensively only by a statistical approach. Our 
sample of stars is not large enough for this. We note however, 
that we do not find a significant trend of carbon abundances 
with evolutionary age for these apparently slow-rotating stars 
(only rSco is suggested to be a real slow rotator, Donati et 
al. 120061) . The abundances derived may not correspond to 
the pristine values nonetheless, but it is unlikely that all the 
stars are very fast rotators seen pole-on. Small corrections by 
up to — i-0.05dex per star (depending on evolution stage) are 
predicted by theory for objects of average rotation. We estimate 
a correction of +0.03 dex to the average value of our sample 
(considering a small increase of the uncertainty), in order to 
derive an average pristine abundance of e(C) - 8.35 ± 0.05. 

The average pristine carbon abundance from early B-type 
stars in the solar neighbourhood is compared with other indi- 
cators in Table|4] The value is not only in agreement with the re- 
vised solar abundance but also with the gas-phase abundance de- 
rived for the Orion nebula. The finding of a uniform abundance 
matches with a homogeneous ISM abundance. A disagreement 
is found only with abundances from young F & G-type stars. 
This may be resolved when hydrodynamic 3D-analyses of late- 
type stars become routine, as the discrepancy is similar to the 
difference between the old and the revised solar standard, which 
is the result of such an improved analysis. The conclusion of 
Sofia & Meyer (2001) that B-type stars are not reliable prox- 
ies for present-day abundances in the solar neighbourhood may 
need a revision in view of the present results. 

Despite the small sample size analysed so far, our accurate 
results indicate that the sub-solar average value and the large 
scatter of C abundances in early-type stars found in previous 
non-LTE studies could be mostly a consequence of systematic 
uncertainties. These may be easily introduced by the choice of 
inappropriate atomic data and/or stellar parameters, as we have 
shown. However, a confirmation of our findings from analysis of 
a larger sample of stars is required, also for other elements. 

7. Summary and Conclusions 

The motivation of this work was the solution of a long-standing 
problem in stellar astrophysics: the reliable determination of 
carbon abundances from early-type stars. For this purpose we 
constructed a sophisticated Cii/iii/iv model atom for non-LTE 
line-formation calculations based on input atomic data that were 
carefully selected by an empirical calibration process. This was 
performed through an extensive iteration scheme that not only 
allowed us to constrain the input atomic data but also simulta- 



neously the atmospheric parameters of our calibration stars, the 
basis for all further studies. The calibration sample consists of 
six bright and apparently slow-rotating early B-type stars in the 
solar neighbourhood, with high-S/N, high-resolution spectra and 
a broad wavelength coverage. 

The self-consistent analysis provides atmospheric parame- 
ters with unprecedented accuracy and with reduced systematic 
error: for effective temperature the uncertainties are as low as 
~1% and for surface gravity ~1Q%. The atmospheric parame- 
ters derived from the carbon ionization equilibria are highly con- 
sistent with our previous spectroscopic analyses: results from a 
non-LTE study of H and He in the visual and in the near-IR, in- 
cluding the establishment of the He i/ii ionization equilibrium in 
the hotter stars (Paper 1), and the reproduction of the spectral 
energy distributions from the UV to the near-IR (Paper 2). 

All C lines considered, from up to three ionization stages, 
indicate similar abundances. The linelist includes 40 transitions 
suitable for analysis over a wide wavelength range. In particu- 
lar, the strongest features, of highest importance for extragalac- 
tic applications, are consistently modelled. The statistical Icr- 
uncertainties of the C abundance in each star are of the order 
~0.05-0.10dex. We estimate the systematic Icr-uncertainties to 
be ~0. 10-0. 15 dex. The average carbon abundance from C ii/iii/iv 
in the sample stars is highly uniform, e(C) = 8.32+0.04. 

These results of unprecedented precision provide important 
constraints for stellar evolution and the chemical evolution of 
the Galaxy, despite the small sample size: i) They suggest that 
carbon depletion due to rotational mixing in the course of stel- 
lar evolution is small for stars without excessive initial angular 
momentum (by <0.05 dex) up to the giant stage, in agreement 
with theoretical predictions, ii) In consequence they suggest that 
the present-day carbon abundance in the solar neighbourhood 
is higher and much more homogeneous than indicated by pre- 
vious work on early-type stars. This is consistent with the uni- 
form abundance found in studies of the ISM and with predictions 
from models of the chemical evolution of the Galaxy. Moreover, 
our re-evaluated stellar value is in agreement with the gas-phase 
abundance derived for the Orion H ii region and the recently re- 
vised solar abundance. 

This work shows that one should not underestimate the im- 
portance of a careful choice of input atomic data for non-LTE 
analyses. Not only accurate data for radiative transitions are re- 
quired but also for collisional transitions. The carbon abundance 
analysis also turned out to react highly sensitively to the choice 
of atmospheric parameters, in particular to effective tempera- 
tures. If not accounted for properly, both factors may result in 
systematic errors in the interpretation of observed spectra, which 
cannot be reduced by statistics, i.e. an increase of the sample size 
of stars to be analysed. Detailed comparisons and calibrations of 
models with observed spectra are of highest importance, as these 
constitute the only empirical constraints for quantitative analy- 
ses. Only then can stellar and Galactochemical evolution models 
be verified in a meaningful way, and the studies be extended to 
stars in other galaxies. 
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Fig. 17. Examples of our line profile fits (black line) to the observed high-quality spectrum (grey line) of t Sco. 
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Fig. 19. As Fig. [181 but for HR 1861. 
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Fig. 20. As Fig.[l8l but for HR 2928. 
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Fig. 21. As Fig.[T8l but for HR 5285. 
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Fig. 22. As Fig.[T8l but for HR 3468. 



